Spin Hall effect in clean two dimensional electron gases with Rashba spin-orbit 

coupling 
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We study the spin polarization induced by a current flow in clean two dimensional electron gases 
with Rashba spin-orbit coupling. This geometric eflect originates from special properties of the 
electron's scattering at the edges of the sample. In wide samples, the spin polarization has it largest 
value at low energies (close to the bottom of the band) and goes to zero at higher energies. In this 
case, the spin polarization is dominated by the presence of evanescent modes which have an explicit 
spin component outside the plane. In quantum wires, on the other hand, the spin polarization 
is dominated by interference effects induced by multiple scattering at the edges. Here, the spin 
polarization is quite sensitive to the value of the Fermi energy, especially close to the point where a 
new channel opens up. We analyzed different geometries and found that the spin polarization can 
be strongly enhanced. 

PACS numbers: 72.25.Dc,73.23.-b,73.23.Ad,71.70.Ej 



I. INTRODUCTION 

The possibility to create and manipulate spin currents 
or to induce spin polarization at the edges of a semi- 
conductor by simply applying an external charge current 
has created great expectation in the field of spintronics — 
a fast developing area aimed to build up new technolo- 
gies based on the manipulation of the electron's spin.^ 
The phenomenon of current induced spin polarization 
(CISP) — also called spin Hall effect (SHE) when referred 
to the perpendicular polarization in two dimensional 
systems — was predicted some time ago by D'yakonov 
and PereliS However, recent theoretical predictions^i^ 
and, more lately, the experimental observation of the 
effec^SiLSiS have renewed the interest on the subject and 
triggered a tremendous among of work during the last 
two years. Although some of the initial controver- 
sies has been settled, some aspects of the problem are 
not completely understood yet. 

Current induced spin polarization may occur in sys- 
tems with strong spin-orbit (SO) coupling, "^^ where the 
electron's momentum and its spin are coupled.**^ The 
simplest setup for observing a CISP consists of a bar- 
shaped sample connected to two terminals. When a 
charge current flows through the system, a spin polar- 
ization develops at the lateral edges of the sample with 
different signs on each side. The mechanisms leading to 
this polarization can be of different nature. In the con- 
text of the SHE, they have been classified as extrinsic 
or intrinsic. The former occurs in dirty semiconductors, 
when the presence of impurities generates a spin depen- 
dent scattering which deflects different spin projections 
in different directions. ^•'^ The latter, on the other hand, 
is due to the presence of an intrinsic SO coupling, which 
is non-local in space and results from the breaking of the 
inversion symmetry of the sample or, in the case of an 
heterostructure, of the confining potential4£ This is the 
mechanism underlying the prediction of dissipansionless 



spin currents in p-doped semiconductors^ and in two- 
dimensional electron gases (2DEG) with Rashba spin- 
orbit couplingi^ 

In the last case, it was argued that the spin Hall con- 
ductivity, defined as the ratio between the transverse spin 
current (generated by the charge current flow) and the 
external electric field, have the universal value e/Sn^ 
Both, the fact of this value being independent of the 
scattering time, and the fact that the spin current is 
not conserved in systems with SO, lead to some contro- 
versy. Furthermore, because of the non-conservation of 
the spin current, the very existence of CISP is not ob- 
vious. Later works, however, have shown that the spin 
Hall conductivity, as obtained in infinite homogeneous 
systems, is very sensitive to disorder and that, in the 
case of Rashba coupling, it vanishes even in the weak 
disorder limit. ^'^i^^i^'^i^^i^'^i^^i'^^i'^" It is important to em- 
phasize at this point that the finite size of the sample was 
not taken into account explicitly in early work but have 
been included in more recent onesi ^^i^'^i^^'^^i?^!?^!?^'?^!?? 

Recent experiments, on the other hand, succeed in ob- 
serving CISP in semiconducting heterostructureaiS and 
in GaAs thin filmsi^ The optically detected magnetiza- 
tion at the edge of a finite sample clearly shows the exis- 
tence of the effect. In addition, CISP at the corner of a 
L-shaped sample was also observed,^ showing the effect 
is not restricted to a straight sample configuration. The 
nature of the observed SHE is, nevertheless, still unclear. 
Arguments in favor of an intrinsic SHE have been put 
forward in Refs. [^35j| while other s ha ve argued in favor 
of an extrinsic explanation, Refs. [dQllsH 

The debate on the physical origin of the CISP triggered 
a number of numerical studies. ^"^■^^■^^i^'^i^^i'^^i^^i^^i^^i^^i^''' 
For the case of clean 2DEG with Rashba SO-coupling, 
it was shown in Ref. that spin polarization near the 
sample edge can develop kinematically: a ballistic current 
that unbalance the occupation of positive and negative 
velocity states leads to a net polarization perpendicular 



to the 2DEG^i — similar results were obtained in Ref. 
for the case of narrow samples. This effect is due to the 
non-trivial structure of the wave functions as result of 
the boundary condition. It is important to emphasize 
that in this clean system, the CISP does not result from 
the presence of a 'bulk' spin current. Therefore, when 
interpreting the CISP in ballistic system as resulting from 
the action of a 'force', it is important to make clear that 
the boundary condition plays a critical lolem^^ In this 
sense, it is worth pointing out that in the case of a bulk 
system, a particle with spin 'up' in the z-direction, will 
not drift but rather follow an oscillatory trajectoryiiSi 

Within this context, it is then important to clarify 
what is the linear response of a ballistic system, what 
is the sample size dependence of the CISP or how the 
geometry can generate different profiles for the spin po- 
larization. Here, we address these questions by studying 
the spin response of nanoscopic systems to an external 
bias using the Keldysh formalism^^i^ 

The paper is organized as follows: the main features 
of the electron's scattering at the surface is revisited in 
Section II for the case of semi-infinite 2DEGs. Finite size 
samples and geometry effects are analyzed in Section III. 
We summarize and conclude in Section IV. 



II. SPIN POLARIZATION IN WIDE SYSTEMS 

A. Semi-infinite systems 

Our starting point is a system with the simplest ge- 
ometry able to show CISP: a semi- infinite 2DEG. Since 
we consider the ballistic regime, the reflection at the 
sample's edge is the only scattering process. As we 
discuss below, this scattering process has very unusual 
propertiesi224SiiLi2 The most relevant are: (i) the ap- 
pearance of evanescent modes localized at the edge and 
(ii) the mixing of bulk states from different bands. As 
we shall see, this geometric effect is enough to lead to 
spin polarization in the presence of an external charge 
current. 

Let us then consider a 2DEG with Rashba SO-coupling 
described by the following Hamiltonian 

H^^^ + Vir) + Hso (1) 
Zm* 

where m* is the effective mass and V{r) is the lateral 
confining potential. The last term accounts for the SO 
intcraction)^2iS 

Hso'^^^iPy^x-Px^y) (2) 

with a the Rashba coupling constant. It is read- 
ily seen from the above equations that the SO cou- 
pling acts as a strong in-plane magnetic field, propor- 
tional to the momentum an perpendicular to it, B-mt = 
{gfiB/2)~^a/h{py,—px). The presence of this effective 



2 



a) K 




FIG. 1: (color online) Schematic representation of the scatter- 
ing process at the edge. Panels a) an b) correspond to incident 
waves with kx <k+. In those cases, two waves are reflected, 
one on each band. Panel c) shows the case of > k+. Here, 
there is only one reflected plane wave but an evanescent mode 
localized at the edge appears. The spin of the localized mode 
has a nonzero component both in the y and the z direction. 

field breaks the spin degeneracy and leads to two non- 
degenerate conduction bands. In an homogeneous 2DEG 
{V{r) = 0) the eigenstates are given by 

with e"^ = {ky — \kx)/k and A the area of the system. 
Note that the electron's spin points in the direction of 
Bint, j-c. it is perpendicular to the momentum (see Fig. 
0. The corresponding eigenvalues are 

e±{\<.) = n^k'^/2m*±ak. (4) 

Then, for a given energy e, there are two characteristic 
wavevectors, and fc_ (with < fc_). We now intro- 
duce the sample's edge by the following potential 

r for,>0 
^ ^ 1 oo for y < ^ ' 

Because of the translational invariance of the system 
along the x-axis, an incident wave (k^,ky),r){f), where 
?7 = ± is the band index, is reflected conserving the x- 
component of the momentum, k^ ■ The continuity of the 
wavefunction at the edge requires ^/'fe^^(a;,0) = 0, where 
'ipk,riif) is the total wavefunction. Since two states on 
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FIG. 2: a) Constant energy surface in the K — kx plane for 
e = and £7^0. © and © indicates spin 'up' and spin 'down' 
in the 2-axis, respectively, b) Schematics of the eigenstates 
of a wide sample for kx > k+. Notice that each states has 
opposite sign of the spin on opposite sides of the sample and 
that the sign of the spin depends on the sign of kx- 
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FIG. 3: Current induced spin polarization in a semi-infinite 
2DEG for different values of Ep. We used a = lOmeVnm, 
m* = 0.068mo and eV — O.lSmeV. The inset correspond to 
a convolution with a Gaussian of rms — 0.25^m. The CISP 
shown in the bottom panels were multiplied by the factors 
indicated in the figure for the propose of comparison. 



the same band with the same have different spin ori- 
entations, the boundary condition can only be satisfied 
if the reflected wave is a linear combination of the two 
modes ^'(fc,,-fc^),,,(r) and ^'(fe, -fe;) -r,(r) as schemati- 
cally shown in Figs. QJi and This superposition of 
plane waves with different spin generates oscillations of 
the spin density. This interference effect have interest- 
ing properties that manifest more clearly in narrow wires 
(see next Section). For the geometry we are considering 
here (or for wide samples), it is more convenient to ana- 
lyze first the case of large incident angles. It is clear from 
FigCl: that for > k+ = [{2^:^6/711,* +a^)i -\a\]m* /h^, 
there are no propagating modes in the band. How- 
ever, for those values of k^, the Schodringer equation 
admittes an evanescent wave solution given by 



*cv(r) = ^e"^-^^e^^( 



(6) 



where 



1,2 1,2 



\a\k^a = ia{kx — K)b , (7) 



C is a normalization constant and |ap + |6p = 1. A 
generic energy surface in the K-kx space is shown in 
Fig|21 While in bulk this solution does not have any 
physical meaning, in the presence of a surface it must 
be included in order to satisfy the boundary condition. 
The total wavefunction is then a linear superposition of 
*(fc,,fe^) _(r), _(r) and -^k^^evir). The evanes- 

cent modes have two interesting properties which are key 



ingredients in leading to CISP in wide samples — the solu- 
tions for smaller incident angles share some of those prop- 
erties, as we discuss below, but are less obvious. First, 
we notice that the evanescent modes have an explicit spin 
component in the z-direction. This is clear from the fact 
that 



\a\\_ kx-K 
|6P kx + K^ 



(8) 



Moreover, the sign of the spin projection depends only on 
the sign of k^ and K . Since the sign of K is fixed by the 
condition that the wavefunction remains finite for y ^ 00 
{K < in our case), the sign of the spin projection is 
giving only by the sign of k^- Then, electrons with fc^, > 
have spin 'up' [a? /h"^ > 1) and those with kx<0 have spin 
'down' {a? /h"^ < 1). Note that for e = the evanescent 
states are fully polarized, a = 1 for k^ >0 and 6=1 for 
kx <0 and that they are the only states available since 
fc+ = 0. It is worth mentioning that in the case of a wide 
but finite sample, the solution with K < dominates 
at the lower edge while the one with K > does it at 
the upper edge. Then, for a given state (value of kx) 
the spin projection have opposite signs in opposites sides 
of the sample. This is schematically shown in Fig. [St. 
Notably, the above results do not depend on the sign of a. 
This result is consistent with the fact that the effective 
'force'iSSiiS acting on a wave package is proportional to 

This simple analysis of the evanescent modes shows 
that there is an intrinsic asymmetry that arises from the 
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combined effect of the breaking of translational invari- 
ance and the intrinsic chirahty introduced by the spin- 
orbit interaction, Vso {p x a) ■ z. This is all what 
is needed to justify the appearance of CISP at e = 0. 
For higher energies, the states without evanescent modes 
must be included. Although each of those states do not 
have a well defined spin component in the z-direction (it 
oscillates as a function of y) , the sum of all the states with 
kx < fc+ leads to a non-zero contribution which has the 
opposite sign that the one corresponding to the evanes- 
cent modes (see Fig0J. Despite it is hard to prove an- 
alytically that such contribution is non-zero, it is rather 
straightforward to show that each state has the property 
that the z-component of the spin density changes sign 
when kx does it. 

In order to calculate the CISP in this ideal geometry, 
we consider a biased 2DEG where the voltage drops at 
the contacts (not included in the calculation) so the elec- 
tric field is zero inside the sample. Then, in the energy 
interval [Ep—eV/2, Ep+eV/2], the carriers injected from 
the left contact, which has the highest chemical poten- 
tial, occupied only the states with kx > Oi^iiS, The CISP 
is then given by 



(9) 



fc,?);fcx>0 



where {az)k.r, = ipl jj{r)a-;,ipk.r,{r) with Fk^,,{EF,eV) = 
e{EF + eV/2-e,^{k))~e{EF-eV/2-e,j{k)). In what 
follows we calculate {Sz{y)) in linear response by numer- 
ical integration of the Schrodinger equation using finite 
differences. The advantage of this method is twofold: a) 
it is not restricted to a square well confining potential; 
b) it can be easily generalized to include the contacts 
(see next section). All quantities presented below are ob- 
tained from the one particle propagators which are cal- 
culated using a continues fraction method (see Ref.js^ 
for details). 

Figure 131 shows {Sz{y)) for different values of Ep and 
for a = lOmeVnm. The insets correspond to the same 
results convoluted with a Gaussian. It is clear from the 
figures that the contribution from the evanescent modes 
is large at low energies and that the CISP goes to zero 
when the energy is increased. The decay of the CISP 
with energy is shown in Fig. 0] where we plotted the 
integrated spin polarization, defined as 



1 = 



(10) 



as a function of Ep. Both, the total integrated spin po- 
larization (Xt, solid dots) and the one coming from the 
evanescent modes (Xov, open dots) are shown. Clearly, 
the decay is well described by a power law. The fact that 
Xt is smaller than X^v indicates that the spin polariza- 
tion of the evanescent modes is opposite in sign to the 
rest of the modes. The decay of lev can be fitted very well 
as lev = h{Ep-\-a'^ra* I2h)~i where 6 is a normalization 
constant 




Fermi energy [meV] 

FIG. 4: Integrated spin polarization T as a function of iJp for 
a = lOmeVnm. The solid dots show the total spin polarization 
while the open dots show the contribution of the evanescent 
modes. Notice X decays with iJp as a power law (the lines are 
a guide to the eyes). 



B. Wide wires 

The case of a wide but finite sample can be analyzed 
in a very similar way. The confinement potential is now 
given by 



V{v)-- 



for < y < 
oo otherwise 



(11) 



where Ly is the width of the wire. The only difference 
here is that a second boundary condition must be im- 
posed, namely tpk,T){x,Ly) — 0. It is clear from Fig. Q] 
that for a given energy, the four states having the same 
kx are mixed. Here again the evanescent modes must 
be included when kx > kj^ (in this case, both solutions 
with K > and K < Q are required) . Figure |S1 shows 
the CISP as a function of the width of the sample. No- 
tice the magnitude of the effect, when averaged over a 
fixed area, is reduced as the sample becomes wider and 
that the CISP is more localized at the sample's edge. 
Besides the oscillating pattern on the scale of the Fermi 
wavelength, there is a longer modulation set by the SO 
coupling, Lgo = 7rfi,^/m*a, which is independent of the 
system's size. 

So far, we have shown that CISP is possible in clean 
systems with SO coupling due to the special properties 
of the electron scattering at the sample's edge. In very 
wide samples, the effect is essentially confined to a region 
of order Lso near the edges and its magnitude decay as 
the Fermi energy increases. As the sample becomes nar- 
rower, however, multiple scattering at the edges leads to 
a different profile for the CISP. It is important to em- 
phasize that there is no conceptual difference between 
both cases and that the picture described so far is the 
basic ph ysic s underlying the mesoscopic SHE described 
in Refs, 



l ysic s underly: 
|2l2l3l53 |. 
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FIG. 5: (color online) Current induced spin polarization 
for diflFerent sample's width Ly = 5,2, 1,0.5/im and for a — 
lOmeVnm, £p = 10meV and ey = 0.13meV. Notice the CISP 
is modulated on a length scale Lso =7rft^/m*a. 



III. SPIN POLARIZATION IN NARROW 
SAMPLES 

Let us now discuss the CISP in finite samples. In order 
to do that, we expHcitly include the leads in our calcu- 
lation. For simplicity, we consider leads with the same 
microscopic parameters than the sample except for the 
SO coupling parameter a which is taken to be zero inside 
the leads. Since the SO parameter a has a spacial varia- 
tion at the sample-lead interfaces, which is assumed to be 
only in the z-direction, the Hamiltonian includes a term 
proportional to —iayda{x)/dxJ^^ As mentioned before, 
to be able to describe systems of arbitrary shape, and in 
order to include the leads, it is convenient to reduce the 
continuum effective mass model to a tight binding model 

H — £(tcJjctC„o- ^ ^nrnCn^C^^ + h.C. 



t 



n 

-A„,„+2 (cl|C(„^2^-cJjj^C(„_,^2)l^ I + h.c.{l2) 

where cj^^ creates an electron at site n with spin az—cr 
and energy = 4t, tn^m = t = /2m*aQ for neighbor- 
ing sites and ao is the lattice parameter. The SO cou- 
plings are defined as follow: \n,n+x = \i,n+y = A = a/2ao 
deep inside the sample, A„,n+2 = (a„-|-a„^„+2)/4ao and 
A„,n+5 =Q!n/2ao at the sample-lead interface and \i^m = ^ 
otherwise. The summation is carried out on a square lat- 
tice where the coordinate of site n is rn — rixX + Uyy with. 
X and y the unit lattice vectors in the x and y direc- 
tions, respectively. Unless otherwise mentioned, we use 
a smooth function to turn on the SO coupling a in order 
to avoid stron g reflections at the sample- lead interfaces. 



a=0 




a=0 


Hi = eV/2 




Hr= -eV/2 



FIG. 6: Schematic representation of a finite quantum wire. 
The SO coupling is non-zero in the dark area. 



When a voltage bias V is applied between the leads, 
an electric current flows through the sample. In this sta- 
tionary but out of equilibrium state, the SO coupling 
generates a pattern of local spin polarization. Within 
the Keldysh formalism, the spin pattern is given hyS^i^ 

(5a(r„)) = ^Tr[aaG<„(t = 0)] (13) 

where the trace is taken on the spin variables, Ca 
with a = x,y, z are the Pauli matrices and the lesser 
Green function (t) is a 2 x 2 matrix with elements 
i(cjjg.(t)c^g., (0)). In what follow we consider small bias 
voltages and linearize Eq. (|13|l . The Fourier transform 
m of the lesser Green function is given byMii^ 

G<„ (e) - i [G^ (a) {hie)r^ + fnie)T^) G'^ (e)]_ 



.eVdfie) 



(14) 



where G'-(e) = [e-i?-S£-S^]-iand G" (e) = [G'^(£)]t are 
the retarded and advanced Green functions, respectively, 
I]£ and are the self-energies due to the left and right 
leads, f{e) is the Fermi function and F'' =i(SJJ— SJ^) with 
r] = L, R . The local spin polarization in the zero temper- 
ature limit and in the absence of an external magnetic 
field, is then given by 

(5a(r„)) = ^Tr [a, {G^iE^)iT^ - r«)G''(i?F)}]„„ 

(15) 

In what follows we use this expression to evaluate the 
CISP in different geometries. 



A. Quantum wires 

We consider a two-terminal wire as shown in Fig. 
For the case of symmetric samples, as the one considerer 
here, all components of the CISP have a well defined 
symmetry. 



{SyiR))v 
{Sz{R))v 



-{S,{-R))-v 
-(Sy{-R))^v 
{SA~R))-v 



(16) 



where R = {x, y) is the spatial coordinate measured from 
the sample's center. Moreover, a 7r-rotation of the system 



100 



200 



FIG. 7: (color online) Current induced spin polarization in 
a quantum wire with Lx = 400ao, Ly = 50ao, -Ef ~ 5meV 
and a = lOmeVnm. (a) Spatial profile of the SO coupling 
parameter a; (b) {Sx{x,y)). Notice it is very small inside the 
sample; (c) {Sy{x,y)); and (d) {Sz{x,y)). Panel (e) shows the 
convolution of {Sz{x,y)) with a Gaussian (we multiplied by a 
factor 5 to keep the same color scale). 



along the x-axis gives 

{Sx{x,y))v,a = {Sx{x,~y))v-a^-{Sx{x,-v))v.a 
{Sy{x,y))v,a = -{Sy{x,-y))v-a^{Sy{x,-y))v.a 

{Sz{x,y))v,a = -{Sz{x,-y))v.-a^~{Sz{x,-y))v.a 

(17) 

where the last equality is due to the symmetry of the 
Hamiltonian, H{a, ax, (Jy) = H(—a, —<Tx, — Cy). Similarly 
a TT-rotation along the y-axis gives 

{Sx{x,y))v,a = ~{Sx{~X,y))-V-a^{Sx{-X,y))-V,a 

{Sy{x,y))v,a = {Sy{-x,y))-v-a = -{Sy{-x,y))-v,a 
{Sz{x,y))v.a = -{Sz{-x,y))^v,~a = -{Sz{~x,y))^v,a 

(18) 

In linear response, these relations imply that;^ 

{Sx{x,y)) = -{Sx{-x,y))=-{Sx{x,-y)) 
{Sy{x,y)) ^ {Sy{~x,y)) = {Sy{x,-y)) 
{Sz{x,y)) = {S,{-x,y))=-{Sz{x,~y)). (19) 

Note that in systems with translational invariance (or in 
the case of very large samples) , the first equation implies 



FIG. 8: (color online) Same as in Fig. |7|but for =500ao, 
Ly = 100ao and i?F=4.8meV. The polarization densities were 
multiplied by the factors indicated in the figure to keep the 
same color scale). 



{Sx{x, y)) =0, in agreement with the results of the previ- 
ous section. Any deviations from the symmetry relations 
shown in Eq. (|19|l , such as a larger spin polarization near 
the drain contact, can only arise from non-linearitiesi^ 

The transverse structure of the different components 
of (S) is determined by three characteristic lengths: the 
Fermi wavelength Xp, the SO length = 7rfi,^/m*a = 
aoTTt/X and the sample width Ly. For wide samples 
{Ly/Lso ^ 1) it has been shown in the previous sec- 
tion that the CISP presents transverse oscillations with 
the two characteristic lengths Xp and Lso- The resulting 
structure shows a beating of the two lengths with a net 
spin polarization at the sample border. 

Results for narrow samples {Ly/Lso ^ 1) are shown 
in Figs. [71and|Hl AH the symmetries mentioned above 
are apparent in the figures. The profile of the Rashba 
parameter a is shown in the top panels. We notice that: 
(i) {Sx) (panel (b)) is zero away from the sample-lead 
interfaces; (ii) (5*1,) is the largest spin component hav- 
ing a defined sign throughout the channel; and (iii) {Sz) 
(panel (d)) oscillates with a dominant sign on each side 
of the wire. To emphasize the later behavior, we made a 
convolution of {Sz) with a Gaussian with an rms of three 
lattice parameters (ISnm). The convoluted spin density 
(panel (e)) clearly shows the spin polarization with oppo- 
site signs at the two sides of the wire. As oppose to the 
wide sample case, however, the transverse profile shows 
strong oscillations whose magnitude is comparable with 
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100 200 300 400 
X W 

FIG. 9; (color online) Same as in Fig. |H1 but for the case 
the SO coupling a is turned on abruptly at the sample-lead 
interfaces. Notice (Sx) inside the sample and that the 
three components of the spin density show oscillations along 
the a;-axis. 



the value of the magnetization close to the edges. 

It is interesting to note that the penetration of the 
x-component inside the sample depends on the sample- 
lead interface. When the Rashba parameter a is turned 
on abruptly, the interference induced by the now strong 
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4.5 5,0 5.5 
E [meV] 

FIG. 10: (color online) a) Local density of states for x — 
260ao; b) Energy dependence of {Sz{x — 260ao,y)). Note 
that the number of oscillations in the j/-axis increases with 
the number of channels; and c) conductance as a function of 
Ep. Parameters as in Fig. |H] 
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FIG. 11: (color online) CISP in a 'T' shaped sample. Elec- 
trons flows from left to right. The left and right lead are 
connected to the system at a; = l and a: = 530, respectively. 



scattering at both interfaces leads to oscillations of (Sx) 
along the x-axis — such oscillating pattern is also visible 
on the other two spin components. This is shown in Fig. 
1^1 Another important feature is that in long samples ( 
Lx/Lso ^ l)i E^nd when a is turned one adiabatically, the 
finite size results away from the interface are the same as 
those obtained in the previous section for infinite long 
systems. 

The CISP is dominated by the channel that contributes 
the most to the density of states (DOS) at the Fermi 
level. This is shown in Fig. IIUI where the transverse 
structure of (Sz) is plotted as a function of Ep for a fixed 
value of X. Clearly the z-component of the CISP is sensi- 
tive to the position of the Fermi energy and the number 
of oscillations depends on the wavelength of the corre- 
sponding transverse mode. In particular, if E-p lies close 
to the bottom of a channel, the magnetization pattern is 
strongly modified: its magnitude increases and its sign 
can change. This behavior suggest that when large volt- 
ages are applied to narrow samples, and therefore a wide 
energy window around E-p is involved, non-linear effects 
might be important due to the contribution of states at 
the bottom of each channel. It also shows the sensitivity 
of the CISP to the geometry while emphasizes the na- 
ture of the effect: it arises from the interference of the 
different spin states which are mixed by the scattering at 
the interface and the non-equilibrium occupancy of the 
momentum states. 
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FIG. 12: (color online) CISP in a 'T' shaped sample when 
Ef is tuned to the value of a resonant state in the cavity. 
Notice the large amplification and localization of the induced 
polarization. 



B. 'T' and 'L' shaped samples 

As mentioned in the previous sections, in the case of 
small samples, both the magnitude and the profile of the 
CISP rehes on the geometry of the system. However, 
so far we have discussed only the wire geometry. In this 
section we will discuss two different geometries: a 'T' and 
a 'L' shaped sample. 

Figures ^2 E^nd ^] show the three components of the 
CISP for the case of a 'T' shaped sample: it corresponds 
to a quantum wire with a open cavity attached to it. 
The cavity is large enough to originate a (resonant) state 
which is substantially localized around the cavity. Figure 
Illl corresDond to an arbitrary value of Ep while in Fig. 1121 
it is tuned to correspond to a resonant state. Notice that 
all components are different from zero. While in both 
cases there is some additional polarization in the vicinity 
of the cavity, the latter case shows a strong enhance- 
ment. We note that because the cavity is at the center 
of the system, some of the symmetry relationships shown 
in Eq. H19|l are still valid. These results suggest that an 
adequate engineering of the sample geometry might al- 
low to design different profiles for the spin polarization, 
such as focusing the spin polarization on a given spot, 
for instance. Notably, the amount of spin polarization 
is very large — in systems without SO a Zeeman field of 
_B ~ 1 — lOT, would be required to achieve a similar 
polarization density. 

Figure [T51 shows the corresponding results for the 'L' 
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FIG. 13: (color online) CISP in a L-shaped system. Notice 
the symmetry between the transverse and longitudinal spin 
components. 



shaped sample. The presence of the 'corner' induces some 
additional spin polarization as compared to the wire ge- 
ometry. The symmetry between (Sx) and (Sy) in the two 
'legs' of the 'L' is apparent from the figure — to emphasize 
this, we have plotted —{Sy) instead of {Sy). Because of 
the strong scattering at the corner, the longitudinal spin 
component does not decay to zero, even far from it. This 
is similar to the case of an abrupt turn on of the Rashba 
parameter. 



IV. SUMMARY AND CONCLUSIONS 

We have showed that current induced spin polariza- 
tion is possible in ballistic 2DEGs. The effect is purely 
geometrical and arises from the electron scattering at the 
sample boundary. For this reason, the CISP is strongly 
dependent on the shape of the sample while some spin 
components can be very sensitive to the sample-lead in- 
terface. In addition, the CISP is quite sensitive to the 
value of the Fermi energy. This is true in both wide and 
narrow samples, though for different reasons. In wide 
samples, the CISP contains two distinct contribution: 
one comes from the interference of the two components 
of the scattered wave, which are two bulk Bloch states 
from different bands, and the other from the evanescent 
modes that appear at the boundary. Both contributions 
have different signs and so tend to cancel each other. At 
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low the evanescent modes clearly dominate so there 
is a large CISP. When the value of Ep is increased, the 
cancellation between the two contributions is better and 
the CISP goes to zero. One interesting feature of the 
CISP for this case is that it does not decay monotoni- 
cally to zero as one goes from the boundary to the bulk 
of the sample but oscillates. Then, even on a given side 
of the sample, the polarization can change its sign. It is 
worth mentioning that because of this strong decay with 
E-p, the phenomena described here cannot account for 
the magnitude of the effect observed in Ref. 0. In such 
device, the presence of a strong disorder and an electric 
field seems to play an important role. One possible expla- 
nation, within our intrinsic model, is that disorder might 
rapidly destroy the phase of the scattered wave without 
affecting the evanescent states too much, therefore unbal- 
ancing the cancellation of the two contributions. In any 
case, we believe, that the special properties of the scat- 
tering at the surface plays a crucial role in originating 
the CISP in those samples. 

The case of narrow systems is a bit different. There, 



the CISP is dominated by the interference of the waves 
scattered at both sides of the sample. Because of the 
strong confinement, the density of states have some struc- 
ture and then, depending on E-p, only a few states might 
dominate the CISP. This is more clearly seen when Ep 
is close to the bottom of a transverse channel or, in the 
presence of a side cavity, near a resonant state. In this 
regime the profile of the CISP very much resembles that 
of a single wavefunction. This opens up the possibility 
to design different profiles for the CISP by modifying the 
geometry of the sample. 
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